Structural and paramagnetic properties of dilute Gai_ x Mn x N 
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Systematic investigations of the structural and magnetic properties of single crystal Gai_ I Mn I N 
films grown by metal organic vapor phase epitaxy are presented. High resolution transmission 
electron microscopy, synchrotron x-ray diffraction, and extended x-ray absorption fine structure 
studies do not reveal any crystallographic phase separation and indicate that Mn occupies Ga- 
substitutional sites in the Mn concentration range up to 1%. The magnetic properties as a function 
of temperature, magnetic field and its orientation with respect to the c-axis of the wurtzite structure 
can be quantitatively described by the paramagnetic theory of an ensemble of non-interacting Mn 3+ 
ions in the relevant crystal field, a conclusion consistent with the x-ray absorption near edge structure 
analysis. A negligible contribution of Mn in the 2+ charge state points to a low concentration of 
residual donors in the studied films. Studies on modulation doped p-type Gai_ a: Mn I N/(Ga,Al)N:Mg 
heterostructures do not reproduce the high temperature robust ferromagnetism reported recently 
for this system. 

PACS numbers: 75.50.Pp, 75.10.Dg, 75.70.Ak, 75.30.Gw, 81.05.Ea 



I. INTRODUCTION 

The search for magnetic semiconductors with a Curie 
temperature Tq above room temperature (RT) is cur- 
rently one of the major challenges in semiconductor 
spintronics. 1-3 In single-phase samples the highest Curie 
temperatures reported are ^190 K for (Ga,Mn)As. 4,5 
The magnetic ordering in these materials is interpreted 
in terms of the p — d Zener model. 1,6 This model as- 
sumes that dilute magnetic semiconductors (DMSs) are 
random alloys, where a fraction of the host cations is 
substitutionally replaced by magnetic ions - hereafter 
with magnetic ions we intend transition metal ions - and 
the indirect magnetic coupling is provided by delocal- 
ized or weakly localized carriers (sp-d exchange interac- 
tions). The authors adopted the Zener approach within 
the virtual-crystal (VCA) and molecular-field (MFA) ap- 
proximations with a proper description of the valence 
band structure in zinc-blende and wurtzite (wz) DMSs. 
The model takes into account the strong spin-orbit and 
the k-p couplings in the valence band as well as the in- 
fluence of strain on the band density of states. This ap- 
proach describes qualitatively, and often quantitatively 
the thermodynamic, micromagnetic, transport, and spec- 
troscopic properties of DMSs with delocalized holes. 3 ' 7 

Experimental data for Gai-^Mn^N reveal an astonish- 
ingly wide spectrum of magnetic properties: some groups 
find high temperature ferromagnetism 8-10 with Tq up to 
940 K, 10 however other detect only a paramagnetic re- 
sponse and their results show that the spin— spin coupling 



is dominated by antiferromagnetic interactions. Gen- 
erally, the origin of the ferromagnetic response in Mn 
doped GaN is not clear and two basic approaches to 
this issue have emerged, namely: i) methods based on 
the mean-field Zener model. 1 - according to this insight, 
in the absence of delocalized or weakly localized holes, 
no ferromagnetism is expected for randomly distributed 
diluted spins. Indeed, recent studies of (Ga,Mn)N in- 
dicate that in samples containing up to 6% of diluted 
Mn, holes are strongly localized and, accordingly, Tc be- 
low 10 K is experimentally revealed. 11,12 Higher values 
of Tc could be obtained providing that efficient methods 
of hole doping will be elaborated for nitride DMSs. Sur- 
prisingly, however, electric field controlled RT ferromag- 
netism has been recently reported in Gai-^Mn^N layers, 
with a Mn content as low as x ps 0.25%. 13 These results 
(Tc > 300 K) cannot be explained in the context of the 
p-d Zener model, where the Curie temperature increases 
linearly with the Mn concentration and for x < 0.5% 
Tc should not exceed 60 K; ii) several studies 14-16 ac- 
knowledge the (likely) presence of secondary phases - 
originating from the low solubility of magnetic ions in 
GaN - as being responsible for the observation of fer- 
romagnetism. It has been found that the aggregation 
of magnetic ions leads either to crystallographic phase 
separation, i.e., to the precipitation of a magnetic com- 
pound, nanoclustcrs of an elemental ferromagnct, or to 
the chemical phase separation into regions with respec- 
tively high and low concentration of magnetic cations, 
formed without distortion of the crystallographic struc- 
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ture. It has been proposed recently that the aggregation 
of magnetic ions can be controlled by varying their va- 
lence (i.e. by tuning the Fermi level). Particularly rele- 
vant in this context are data for (Zn,Cr)Te, 17 (Ga,Fe)N, 18 
and also (Ga,Mn)N, 17,19 ' 20 where a strict correlation be- 
tween codoping, magnetic properties, and magnetic ion 
distribution has been put into evidence. 

There is generally a close relation between the ion 
arrangement and the magnetic response of a magneti- 
cally doped semiconductor. Specifically, depending on 
different preparation techniques and parameters, coher- 
ently embedded magnetic nanocrystals [like wz-MnN in 
GaN (Rcfs. 21 and 22)] or precipitates [like e.g. MnGa 
or Mn4N] might in fact give the major contribution to 
the total magnetic moment of the investigated samples. 
In particular, randomly distributed localized spins may 
account for the paramagnetic component of the magneti- 
zation, whereas regions with a high local density of mag- 
netic cations are presumably responsible for ferromag- 
netic features. 23 In the case of low concentrations of the 
magnetic impurity, it is often exceedingly challenging to 
categorically identify the origin of the ferromagnetic sig- 
natures. 

Up to very recently, in most of the reports the obser- 
vation of ferromagnctism or fcrromagnctic-likc behavior 
with apparent Curie temperatures near or above RT, has 
been discussed primarily or even solely based on magnetic 
hysteresis measurements. However, indirect means like 
superconducting quantum interference device (SQUID) 
magnetometry measurements or even the presence of the 
anomalous or extraordinary Hall effect, may be not suf- 
ficient for a conclusive statement and to verify a single- 
phase system. Therefore, a careful and thorough charac- 
terization of the systems at the nanoscale is required. 
This can only be achieved through a precise correla- 
tion of the measured magnetic properties with advanced 
material characterization methods, like e.g. synchrotron 
x-ray diffraction (SXRD), synchrotron based extended 
x-ray absorption fine structure (EXAFS) and advanced 
element-specific microscopy techniques, suitable for the 
detection of a crystallographic and/or chemical phase 
separation. 

The present work is devoted to a comprehensive study 
of the Gai_ x Mn x N (x < 1%) fabricated by metalor- 
ganic vapor phase epitaxy (MOVPE), which was also em- 
ployed by other authors. 13,19 . A careful on-line control 
of the growth process is carried out, which is followed 
by an extended investigation of the structural, optical, 
and magnetic properties in order to shed new light onto 
the mechanisms responsible for the magnetic response of 
the considered system. Particular attention is devoted 
to avoid the contamination of the SQUID magnetome- 
try signal with spurious effects and, thus, to the reliable 
determination of the magnetic properties. Experimental 
procedures involving SXRD, high resolution transmission 
electron microscopy (HRTEM), EXAFS and x-ray ab- 
sorption near-edge spectroscopy (XANES) are employed 
to probe the possible presence of secondary phases, pre- 



cipitates or nanoclusters, as well as the chemical phase 
separation. Moreover, we extensively analyze the prop- 
erties of single magnetic-impurity states in the nitride 
host. The understanding of this limit is crucial when 
considering the most recent suggestions for the con- 
trolled incorporation of the magnetic ions and conse- 
quently of the magnetic response through Fermi level 
engineering. By combining the different complementary 
characterization techniques we establish that randomly 
distributed Mn ions with a concentration x < 1% gen- 
erate a paramagnetic response down to at least 2 K in 
Gai-^Mn^N. In view of our findings, the room temper- 
ature ferromagnetism observed in this Mn concentration 
range 13,19 ' 20 ' 24,25 has to be assigned to a non-random dis- 
tribution of transition metal impurities in GaN. We em- 
phasize that in all reported works on (Ga,Mn)N fabri- 
cated by MOVPE the Mn concentration was well below 
5%. 

The paper is organized as follows: in the next sec- 
tion we give a summary of the fabrication details, in situ 
monitoring of the employed MOVPE process and an 
abridged overview of the characterization techniques, to- 
gether with a table listing the principal properties and 
parameters characterizing the (Ga,Mn)N-based samples 
considered. In Sec. IV the results of the structural anal- 
ysis of the layers by SXRD, HRTEM, and EXAFS are 
reported. These measurements prove a uniform distribu- 
tion of the Mn ions in the Ga sublattice of GaN. Sec- 
tion V is devoted to the determination of the Mn con- 
centration and of the charge and electronic state of the 
magnetic ions. In section VI we give the experimen- 
tal magnetization characteristics of the system obtained 
from SQUID measurements, and interpret the data based 
on the group theoretical model for Mn 3+ ions taking into 
account the trigonal crystal field, the Jahn- Teller dis- 
tortion and the spin-orbit coupling. Finally, conclusions 
and outlook stemming from our work are summarized in 
Sec. VII. 



II. GROWTH PROCEDURE 

The wz-(Ga,Mn)N epilayers here considered are fabri- 
cated by MOVPE in an AIXTRON 200 RF horizontal re- 
actor. All structures have been deposited on c-plane sap- 
phire substrates with TMGa (trimcthylgallium) , NH3, 
and MeCp2Mn (bis-methylciclopentadienyl-manganese) 
as precursors for, respectively, Ga, N and Mn, and with 
H2 as carrier gas. The growth process has been carried 
out according to a well established procedure 26 consisting 
of: substrate nitridation, low temperature (540°C) depo- 
sition of a GaN nuclcation layer (NL), annealing of the 
NL under NH 3 until recrystallization and the growth of 
a ~1 /im thick device-quality GaN buffer at 1030°C. On 
top of these structures, Mn doped GaN layers (200-700 
nm) at 850°C, at constant TMGa and different — over the 
samples series — MeCp2Mn flow-rates ranging from 25 to 
490 seem (standard cubic centimeters per minute) have 



3 



TABLE I: Data related to the investigated Gai-zMn^N. The 
following values are listed: the MeCp2Mn flow rate employed 
to grow the Mn-doped layers, the FWHM of the (0002) reflex 
from GaN determined by ex situ HRXRD, the Mn 3+ concen- 
tration as obtained from magnetization data, the total Mn 
content from SIMS measurements and the thickness of each 
(Ga,Mn)N layer. Letters A and B denote the two different 
growth series 

McMnCp2 thickness of Mn 3 ^ cone. Mn cone, 

sample flow rate (Ga,Mn)N FWHM SQUID SIMS 

number [seem] [ mn] [arcscc] [10 20 cm~ 3 ] [10 20 cm~ 3 ] 
000B 470 <0.06 

025A 25 450 242 0.28 0.3 

050A 50 400 267 0.8 0.6 



0.5 
0.7 

1.1 

1.3 



1.9 



2.4 



been grown. The nominal Mn content in subsequently 
grown samples has been alternatively switched from low 
to high and, vice versa, to minimize long term memory 
effects due to the presence of residual Mn in the reac- 
tor. During the whole growth process the samples have 
been continuously rotated in order to promote the depo- 
sition homogeneity and in situ and on line ellipsometry is 
employed for the real time control over the entire fabrica- 
tion process. The p-type superlattices have been grown 
according to the optimized procedure already reported. 27 
Our MOVPE system is equipped with an in situ Isa 
Jobin Yvon ellipsometer that allows both spectroscopic 
(variation of the optical parameters as a function of the 
radiation wavelength) and kinetic (ellipsometric angles 
vs. time) measurements 28,29 in the energy range 1.5 - 
5.5 eV.In Table I the considered (Ga,Mn)N samples are 
listed together with their specific parameters. 



III. EXPERIMENTAL TECHNIQUES 

A. HRTEM experimental 

HRTEM studies have been carried out on cross- 
sectional samples prepared by standard mechanical pol- 
ishing followed by Ar + ion milling, under a 4° angle at 
4 kV for less than 2 h. The ion polishing has been per- 
formed in a Gatan 691 PIPS system. The specimens 
were investigated using a JEOL 2011 Fast TEM micro- 



scope operated at 200 kV equipped with a Gatan CCD 
camera. The set-up is capable of an ultimate point-to- 
point resolution of 0.19 nm, with the possibility to image 
lattice fringes with a 0.14 nm resolution. The chemical 
analysis has been accomplished with an Oxford Inca en- 
ergy dispersive x-ray spectroscopy (EDS) system. 

B. HRXRD and SXRD experimental 

High-resolution x-ray diffraction (HRXRD) rocking 
curves are routinely acquired on each sample with a 
Philips XRD HR1 vertical diffractometer with a CuK Q 
x-ray source working at a wavelength of 0.15406 nm (~ 
8 keV). A monochromator with a Ge(440) crystal config- 
uration is used to collimate the beam, that is diffracted 
and collected by a Xe-gas detector. Angular (w) and 
radial uj/29 scans have been collected along the growth 
direction for the (002) GaN reflex, in order to gain in- 
formation on the crystal quality of the samples from the 
full width at half maximum (FWHM) of the diffraction 
peak. 

Though being aware that, if great care is exercised, also 
conventional XRD may allow to detect small embedded 
clusters (like in the reported case of Co in ZnO) 30 ~ 32 we 
performed SXRD measurements that gave us the possi- 
bility to additionally carry out in situ annealing exper- 
iments. The experiments have been carried out at the 
beamline BM20 (Rossendorf Beam Line) of the Euro- 
pean Synchrotron Radiation Facility (ESRF) in Greno- 
ble - France. Radial coplanar scans in the 29 range from 
20° to 60° were acquired at a photon energy of 10 keV. 
The beamline is equipped with a double-crystal Si(lll) 
monochromator with two collimating/focusing mirrors 
(Si and Pt-coating) for rejection of higher harmonics, al- 
lowing an acquisition energy range from 6 to 33 keV. The 
measurements are performed using a heavy-duty 6-circle 
Hubcr diffractometer, that is the system is suitable for 
(heavy) user-specific environments (e.g. in our case a Be- 
dome for the annealing experiments was required). 

C. EXAFS and XANES experimental 

The x-ray absorption fine structure (XAFS) measure- 
ments at the Mn-K edge (6539 eV) have been performed 
at the GILDA Italian collaborating research group beam- 
line (BM08) of the ESRF in Grenoble. 33 The monochro- 
mator is equipped with a pair of Si(311) crystals and 
run in dynamical focusing mode. 34 Harmonics rejection 
is achieved through a pair of Pd-coated mirrors with an 
estimated cutoff of 18 keV. Data are collected in the fluo- 
rescence mode using a 13-element hyper pure Ge detector 
and normalized by measuring the incident beam with an 
ion chamber filled with nitrogen gas. In order to mini- 
mize the effects of coherent scattering from the substrate, 
the samples are mounted on a dedicated sample holder 
for grazing-incidence geometry; 35 measurements are car- 
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ried out at room temperature with an incidence angle of 
1° and with the polarization vector parallel to the sample 
surface (E J_ c). For each sample the integration time for 
each energy point and the number of acquired spectra are 
chosen in order to collect ps 10 6 counts on the final av- 
eraged spectrum. Bragg diffraction peaks are eliminated 
by selecting the elements of the fluorescence detector or 
by manually de-glitching the affected spectra. In addi- 
tion, before and after each measurement a metallic Mn 
reference foil is measured in transmission mode to check 
the stability of the energy scale and to provide an accu- 
rate calibration. Considering the present optics setup an 
energy resolution of pa 0.2 eV is obtained at 6539 eV. 

In our context, the EXAFS signal x(fc) is extracted 
from the absorption raw data, n(E), with the viper 
program 36 employing a smoothing spline algorithm and 
choosing the energy edge value (Eo) at the half height of 
the energy step (Sec. VC). The quantitative analysis is 
carried out with the ifeffit/ Artemis programs 37,38 in 
the frame of the atomic model described below. Theoret- 
ical EXAFS signals are computed with the feff8 code 39 
using muffin tin potentials and the Hedin-Lunqvist ap- 
proximation for their energy-dependent part. In order 
to reduce the correlations between variables, the mini- 
mum set of free fitting parameters used in the analysis 
is: AE (correction to the energy edge), Sq (amplitude re- 
duction factor), ARo, ARi (lattice expansion factors, re- 
spectively, for the first Mn-N coordination shell distances 
and all other upper distances) and of Debye- Waller fac- 
tor for the i th coordination shell around the absorber plus 
a correlated Debye model 40 for multiple scattering paths 
with a fitted Debye temperature of 470(50) K. 41 
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FIG. 1: (Color online) Magnetic response at 2, 15, and 200 K 
of a typical (5 x 5 x 0.3) mm 3 sapphire substrate measured for 
both in-plane (darker shade squares) and out-of-plane (lighter 
shade circles) configurations after the application of a cor- 
rection linear in the magnetic field and proportional to the 
magnetic susceptibility of sapphire at 200 K. The magnetic 
moment obtained in this way at 2 K reaches a value that 
would give a magnetization of ~ 1 emu/cm 3 for (a typical) 
200 nm thick layer. Inset: m(H) at 2 K for the same sap- 
phire sample, but without correction. The axes labels are 
the same as in the main panel. This ferromagnetic-like signal 
is isotropic, decreases with temperature and vanishes above 
15 K. 



D. SQUID experimental 

The magnetic properties have been investigated in a 
Quantum Design MPMS XL 5 SQUID magnetometer be- 
tween 1.85 and 400 K and up to 5 T. For magnetic studies 
the samples are typically cut into (5x5) mm 2 pieces, and 
both in- and out-of-plane orientations are probed. The 
(Ga,Mn)N layers are grown on 330 /im thick sapphire 
substrates, so that the TM-doped overlayers constitute 
only a tiny fraction of the volume investigated, and due to 
the substantial magnetic dilution their magnetic moment 
is very small to small when compared to the diamagnctic 
signal of the substrate. Therefore, a simple subtraction of 
a diamagnetic component originating from the sapphire 
substrate and linear with the field only exposes the result- 
ing data to various artifacts related to the SQUID system 
and to arrangement of the measurements, as already dis- 
cussed in Rcfs. 28, 42, and 43. In order to circumvent 
this issue, the magnetic data presented in this paper are 
obtained after subtracting the magnetic response of a 
sapphire substrate with dimensions equivalent to those 
of the investigated sample, independently measured on 
the same holders and according to the same experimental 
procedure. This method, in particular, eliminates a spu- 



rious magnetic contribution that is due to the sapphire 
substrate and is not-linear with the field and, moreover, 
depends on the temperature. As exemplified in Fig. 1 
this extra m(T, H) constitutes a nontrivial and quite siz- 
able contribution to the signal of interest. Additionally 
- as shown in the inset to this figure - the sapphire it- 
self may convey a ferromagnetic response to the signal at 
the lowest temperatures. We have also made sure that 
this method is adequate to eliminate another weak and 
ferromagnetic-like contribution appearing in the data af- 
ter subtracting only the compensation linear in the field. 
This fault is caused by an inaccuracy in the value of the 
magnetic field as reported by the SQUID system, which 
assumes that the field acting on the sample is strictly 
proportional to the current sent to the superconducting 
coil, and disregards the magnet remanence due to the 
flux pinning inside the superconducting windings. 44 This 
remanence in our 5 T system is as high as -15 Oe after 
the field has been risen to H > +1 T and results in a zero 
field magnetic moment of +2 x 10~ 7 emu for our typical 
sapphire substrate. Although the value is small, it lin- 
early scales with the mass of the substrate and it exceeds 
the magnitude of the signal expected from a submicrom- 
eter thin layer of a DMS film. 
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E. SIMS experimental 

The overall Mn concentration in the epilayers has been 
evaluated via secondary-ion mass spectroscopy (SIMS). 
The SIMS analysis is performed by employing Caesium 
ions as the primary beam with the energy set to 5.5 keV 
and the beam current kept at 150 nA. The raster size is 
150 x 150 /urn 2 and the secondary ions are collected from 
a central region of 60 fim in a diameter. The Mn concen- 
tration is derived from MnCs+ species, and the matrix 
signal NCs+ was taken as reference. Mn implanted GaN 
is used as a calibration standard. 



IV. STRUCTURAL PROPERTIES 



loss during annealing. Several radial scans have been ac- 
quired upon increasing the sample temperature in subse- 
quent 100°C steps, and realignment was performed after 
reaching each temperature. The diffraction curves upon 
annealing are shown in Fig. 3, and no additional diffrac- 
tion peaks related to the formation of secondary phases 
have been detected over the whole process. This leads us 
to conclude that the considered (Ga,Mn)N grown with 
Mn concentration below the solubility limit at the given 
deposition conditions is stable in the dilute phase upon 
annealing over a considerable thermal range. This behav- 
ior is to be contrasted with the one reported for dilute 
(Ga,Mn)As, where annealing at elevated temperatures 
provokes the formation of either hexagonal or zinc-blende 
MnAs nanocrystals. 48 ' 49 



As already underlined, it is necessary to ascertain if 
the investigated material contains any secondary phases, 
nanoclusters or precipitates. In this context, it has been 
realized recently 18 ' 21,28,45 that the limited solubility of 
transition metals in semiconductors can lead to a chem- 
ical decomposition of the alloy, i.e. the formation of re- 
gions with the same crystal structure of the semiconduc- 
tor host, but with respectively high and low concentra- 
tion of magnetic constituents. In this work the structural 
properties of the system are analyzed by SXRD, HRTEM 
and EXAFS. 



A. HRXRD and SXRD - results 

In order to verify the homogeneity of the grown 
(Ga,Mn)N layers, conventional XRD measurements have 
been routinely performed. From rocking curves around 
the GaN (002) diffraction peak, the crystal quality from 
the FWHM is verified and we obtain values in the range 
240 to 290 arcsec, indicating a high degree of crystal per- 
fection of the layers. For the (Ga,Mn)N (002) diffraction 
peak we observe a shift to lower angles with increasing 
Mn concentration in the acquired u>/26 scans. This shift 
points to an increment in the c-lattice constant, as it 
has been also reported by others. 46,47 Apart from the 
diffraction peak shift, no evidence for second phases is 
observed in the XRD measurements. These results have 
been confirmed by the SXRD diffraction spectra reported 
in Fig. 2(a), where no crystallographic phase separation 
is detected over a broad range of Mn concentrations. 

The lattice parameters arc determined by averaging 
the values for the two symmetric SXRD diffractions 
(004) and (006) for the c-parameter, and one asymmet- 
ric diffraction (104) for the a-parameter. The variation 
of the lattice parameters with increasing incorporation of 
Mn is presented in Fig. 2b. 

To obtain further information on the solubility of Mn in 
our (Ga,Mn)N layers, in situ annealing experiments have 
been carried out at the ESRF BM20 beamline. Sample 
400 A was annealed up to 900° C in N-rich atmosphere 
at a pressure of 200 mbar, to compensate the nitrogen 



B. HRTEM results 

HRTEM has been carried out on all (Ga,Mn)N layers 
under consideration and independently of the Mn con- 
centration no evidence of crystallographic phase separa- 
tion could be found. This is also confirmed by selected 
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FIG. 2: (Color online) a) SXRD spectra for (Ga,Mn)N sam- 
ples showing no presence of secondary phases over a broad 
range of concentration of the magnetic ions, b) Lattice pa- 
rameters vs. Mn concentration. Values for an undoped GaN 
layer are added for reference 



area electron diffraction (SAED) patterns (not shown) 
recorded on different areas of each sample, where no 
satellite diffraction spot apart from wurtzite GaN are 
detected. In Fig. 4, an example of the HRTEM im- 
ages acquired along the [1010] (a) and [1120] (b) zone 
axis, respectively is given. Through measurements pre- 
viously reported and carried out with the same micro- 
scope, we have been able to discriminate in (Ga,Fc)N 
different phases of Fe-rich nanocrystals as small as 3 nm 
in diameter, and also to detect mass contrast indicating 
the local aggregation of Fe-ions. 18 ' 50 The HRTEM im- 
ages in Fig. 4, in contrast to the case of phase separated 
(Ga,Fe)N, strongly suggest that the (Ga,Mn)N film here 
studied are in the dilute state 

The EDS spectra collected on the (Ga,Mn)N layers 
provide significant signatures of the presence of Mn, as 
evidenced in Fig. 5. The EDS detector and the soft- 
ware we used here identify the Mn elements automat- 
ically, and are sensitive to Mn concentrations as low as 
0.1% (atomic%). The Mn concentration for sample 300A 
and reported in Fig. 5 is found to be 0.18% (atomic%). 



C. EXAFS results 
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FIG. 4: HRTEM images: along [1010] (a) and along the [1120] 
zone axis (b). 
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EXAFS (Rcf. 51) is a well established tool in the study 
of semiconductor heterostructures and nanostructures 52 
and has proven its power as a chemically sensitive local 
probe for the site identification and valence state of Mn 
and Fc dopants in GaN DMS. 53 ~ 57 The crystallinity of 
the films and the optimal signal to noise ratio of the col- 
lected spectra are demonstrated by the large number of 
atomic shells visible and reproducible by the fits below 
8 A in the Fourier-transformed spectra reported in Fig. 6 
for the two representative samples 100A and 490A, rc- 
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FIG. 3: (Color online) In situ SXRD spectra upon annealing 
of sample 400A at different temperatures: no formation of 
secondary phases is detected up to an annealing temperature 
of 900 °C. 



FIG. 5: EDS spectrum of sample 300A, with the identification 
of the Mn peaks [L a (0.636 keV), K a (5.895 keV) and K p (6.492 
keV) ]. 



spectively. In addition, the homogeneous Mn incorpora- 
tion along the layer thickness is tested by measuring the 
Mn fluorescence yield (at a fixed energy of 6700 eV) as a 
function of the incidence angle (not shown). The EXAFS 
response of these two samples is qualitatively equivalent, 
as evidenced in Fig. 6, and this is confirmed by the quan- 
titative analysis. The best fits are obtained by employing 
a substitutional model of one Mn at a Ga site (Mn Ga ) in 
a wurtzite GaN crystal (using the lattice parameters pre- 
viously found by SXRD). The possible presence of addi- 
tional phases in the sample as octahedral or tetrahedral 
interstitials (Mnf, Mnj r ) in GaN or Mn 3 GaN clusters 58 
has been checked by carrying out fits with a two phases 
model. The fraction of the Mn Ga is found to be 98(4) % 
for the pair Mn Ga -Mn 3 GaN, 99(3)% for the pair Mn Ga - 
Mnf and 97(3)% for the pair MnQa-Mn^, respectively. 
With these results we can safely rule out the occurrence 
of phases other than Mn Ga , at least above 5% level. 

The local structure parameters found for the mea- 
sured samples are equivalent within the error bars (re- 
ported on the last digit within parentheses) and aver- 
aged values are given for simplicity. The value of the 
amplitude reduction factor S 2 , = 0.95(5) demonstrates 
the good agreement with the theoretical coordination 
numbers for Mn Ga (considering the in-plane polarization) 
and the correction to the energy edge AEq = —7(1) eV 
supports the XANES analysis (Sec. VC). With respect 
to the lattice parameters previously found by SXRD, 
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FIG. 6: (Color online) fc 2 -weighted EXAFS signal, (a), for 
samples 100A (circles) and 490A (diamonds) with relative 
best fits (solid line) in the region [Rmin-Rmax] and, (b), ampli- 
tude of the Fourier transforms (FT) carried out in the range 
[kmin-kmax] by an Hanning window (slope parameter dk=l); 
the vertical lines indicate the position of the main peaks in 
the FT of the MnJ , MnP and Mn 3 GaN additional structures. 
Their possible presence would be proptly detected as they fall 
in a region free from other peaks. feff8 simulations, (c), for 
the tested theoretical models as described in the text. 



the long range distortion fits within the error (Ai?i = 
0.1(2) %), while the Mn-N nearest neighbors have a 
AR = 2.5(5) % (expansion to 1.99(1) A), in line with 
previously reported experimental results 54-56 and recent 
ab initio calculations. 59 Finally, all the evaluated of at- 
test around the average value of 8(2)-lCT 3 A -2 , confirm- 
ing the high crystallinity of the layers. 



V. PROPERTIES OF HOMOGENEOUS 
SINGLE-PHASE (Ga,Mn)N 

Thus, SXRD, HRTEM and XAFS experiments have 
confirmed the wurtzite structure of the samples, the ab- 
sence of secondary phases, and the location of Mn in 
the Ga sublattice of the wurtzite GaN crystal. Further- 



more, the samples have been investigated to determine 
the actual Mn concentration and the charge state of the 
magnetic ions. 



A. Determination of the Mn concentration 

The depth profiling capabilities of SIMS provide not 
only an accurate analysis of the (Ga,Mn)N layers com- 
position, but allow also to monitor the changes in compo- 
sition along the sample depth. The SIMS depth profiles 
reported in Figs. 7(a) and (b) give evidence that the dis- 
tribution of the Mn concentration riMn in the investigated 
films is essentially uniform over the doped layers, inde- 
pendent of the magnetic ions content as well as that the 
interface between the (Ga,Mn)N overlayer and the GaN 
buffer layer is sharp. This is confirmed by EDS studies, 
which with the sensitivity around 0.1% at. do not provide 
any evidence for Mn diffusion into the buffer. The deter- 
mined total Mn concentration increases with increasing 
McCp2Mn flow rate and the corresponding nyin values 
for the considered samples can be found in Table I. 




Depth [ urn ] 

FIG. 7: (Color online) SIMS depth profiles of Mn, C, O and 
H for the samples: a) 150A and b) 375A. 



B. Energy levels introduced by Mn impurities 

The character of the paramagnetic response of DMS 
depends crucially on the magnetic ion configuration. In 
III-V semiconductors, Mn in the impurity limit substi- 
tutes the cation site giving three electrons to the crystal 
bond. Depending on the compensation ratio, Mn can 
exist in three different charge states and electron con- 
figurations, namely: i) ionized acceptor Mn 2+ , with five 
electrons localized in the Mn d shell. The electronic con- 
figuration of Mn 2+ is d 5 , and the ground level of the 
ion at zero magnetic field is a degenerate multiplet with 
vanishing orbital momentum (L = 0, S = 5/2). The 
magnetic moment of the ion results solely from the spin, 
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and its magnetic contribution can be described by a stan- 
dard Brillouin function for any orientation of the mag- 
netic field. The neutral configuration of Mn 3+ (S = 2, 
L = 2) can be realized in two ways: ii) by substitutional 
manganese d 4 with four electrons tightly bound in the 
Mn d shell; iii) Mn 2+ + hole (d 5 + hole) with five elec- 
trons in the Mn d shell and a bound hole localized on 
neighboring anions. 



C. XANES results 

The XANES spectra allow to determine the redox- 
state of the probed species and give information on the 
structure of the surroundings of the absorbing atom. 60 
Basically, the near edge region resembles the density of 
those empty states, that are accessible via optical tran- 
sitions from the Mn Is shell. 

The goal of our XANES analysis is to determine the 
valence state of Mn and to confirm the MnQ a incorpora- 
tion, in comparison to the findings and analysis carried 
out previously for molecular beam epitaxy (MBE)-grown 
(Ga,Mn)N, and interpreted in terms of Mn 3+ (Refs. 11 
and 61) or Mn 2+ (Rcf. 62). In order to assign the Mn 
valence state, first of all we proceed with a comparison 
of the position of the absorption Mn K-edge to refer- 
ence compounds, like Mn-bascd oxides since we do not 
have available data on Mn-nitridcs. This procedure was 
already adopted by other groups 55 ' 62 but its reliability 
could be questionable; to clarify this point ab initio cal- 
culations are also performed. 

As shown in Fig. 8, the XANES spectra determined 
for two samples differing in Mn concentrations (100A 
and 490A) are identical, confirming a conclusion from 
the SQUID data on the independence of the Mn charge 
state of the Mn concentration. In Fig. 8(a) three spectra 
collected in transmission mode from commercial powders 
of MnO, Mn 2 3 and Mn0 2 , with Mn-valence states 2+, 
3+, 4+, respectively, are used as reference. As seen, with 
the increasing charge state, the edge moves to a higher 
energy, as the accumulated positive charge shifts down- 
wards in energy more the Is Mn shell than the valence 
states, in agreement with the Haldanc- Anderson rule. 

Usually, the edge position is taken at the first inflection 
point of the plot, but in the present case (since the oxide 
spectra exhibit a broad peak that modifies the slope at 
the edge) a better estimate of the edge position is ob- 
tained by considering the energy of the half step-height 
of the background function. In both investigated samples 
this lies at 6550.0(5) eV. For the oxides, their half-height 
energies are determined to be 6545.7(5) eV, 6550.3(5) eV, 
6553.3(5) eV for MnO, Mn 2 3 and Mn0 2 , respectively. 
This would strongly suggest that we deal with Mn 3+ , in 
line with the SQUID results (Sec. VIA). On the other 
hand, taking the position of the inflection points, the 
determined charge state would be 2+, as reported in 
Ref. 62. This demonstrates that, in this case, relying 
only on the edge position to determine the valence state 



is prone to error and strongly depends on the local sur- 
rounding of the probed species. 63 
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FIG. 8: (Color online) Normalized XANES spectra of the 
samples 100A and 490A (points) compared with: (a) the ref- 
erence manganese oxides (MnO, Mn203, Mn02) - the chosen 
edge positions are highlighted by vertical lines; (b) ab ini- 
tio absorption spectra (without convolution) for MnGa in the 
3d 4 and 3d 5 electronic configurations. The inset (c) shows the 
method used to extract the results of Table II, focusing the 
near-edge region for sample 490A with the baseline (Bkg) , the 
relative fit and its components (A1-A4). 

To clarify this point, we look at the pre-edge peaks of 
the XANES lines (Fig. 8(b), (c)). In both probed sam- 
ples, there are two defined peaks below the absorption 
edge, which we label Ai and A 2 , while the edge itself 
shows two shoulders, A3 and A4. In Table II the results 
of Gaussian fits performed by using an arctan-function as 
baseline, are reported. Similar findings were previously 
interpreted 61,64 as indicative of the Mn 3+ charge state. 
The peaks Ai and A 2 correspond to the transitions to Mn 
3d-4p hybrid states, while A 3 and A4 end in the GaN 
higher conduction bands at positions with a high den- 
sity of Ap states. Due to the tetrahedral environment, 
the Mn 3d-levels split in two nearly degenerate e- and 
three nearly degenerate i 2 -levels for each spin-direction. 
The actual position of those states with respect to the 
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TABLE II: Position P, integrated intensity / and full width at 
half maximum W of the Gaussians fitted to the peaks before 
and at the absorption edge. The background function is used 
to normalize the spectra. 







100A 






490A 






P (eV) 
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P (eV) 
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W (eV) 




± 0.2 


± 0.05 


± 0.1 


± 0.2 


± 0.05 


± 0.1 




6538.9 


0.18 


1.5 


6538.9 


0.18 


1.3 


A 2 


6540.8 


0.32 


1.6 


6540.8 


0.33 


1.5 


A 3 


6545.8 


0.70 


3.6 


6545.7 


0.64 


3.3 


A 4 


6548.7 


0.28 


2.1 


6548.4 


0.22 


2.0 



GaN band structure is still a matter of debate, but from 
absorption 65 ' 66 and photolumincscence 67 measurements 
it is known that for the majority spin carriers in Mn 3+ , 
the e- levels lie around 1.4 eV below the ^-levels of Mn 
incorporated substitutionally in GaN, and the i2-level, 
i.e., the Mn 3+ /Mn 2+ state is located about 1.8 eV above 
the valence band. An interpretation of simulations ap- 
plied to x-ray absorption spectra is given in Refs. 61 and 
68, and states that, due to crystal field effects, the 3d- 
and 4p-states can hybridize, making transitions from the 
ls-level to the i2-levels dipole allowed, while the interac- 
tion of the e-levels with the Ap orbitals is much weaker 
and cannot be seen in K-edge XANES. 

In view of the above discussion we explain the physical 
mechanism beyond the observed data considering possi- 
ble final states of the transitions from the Is Mn shell. 
The final state corresponding to the Ai peak is Mn 2+ , 
i.e. a 6 Ai state ( 6 S for the spherical symmetry), consist- 
ing of e 2 ^ and one electron levels. The A2 peak can be 
interpreted as a crystal field multiplct derived from the 
4 G state consisting of e 2 ^ an( ^ lying about 2.5 eV 
higher than the Ai state. Apart from what reported in 
literature, a reason why Ai and A2 are assigned to local- 
ized Mn-states is that from the previous EXAFS anal- 
ysis (Sec. IV C) we obtain an absorption edge value of 
6543(1) eV, between the energies of the A2 and A3 peaks, 
meaning that electrons excited to Aj and A 2 can not 
backscatter at the surrounding atoms, and they are thus 
localized. This assignment gives a valuable information, 
namely, that there is an empty state in the majority-spin 
t2-level confirming that most of the incorporated Mn- 
ions are really in the 3+ valence state, in agreement with 
the conclusions of Refs. 11 and 61. The model explains 
also the presence of only one pre-edge peak in the case 
of (Ga,Mn)As and p-(Zn,Mn)Tc. 61 In those systems we 
deal with Mn 2+ and delocalized holes, so that the final 
state of the relevant transitions corresponds to the Mn d 6 
level, involving only one spin orientation. On the other 
hand, the XANES data do not provide information on 
the radius of the hole localization in (Ga,Mn)N, in other 
words, whether the Mn 3+ configuration corresponds to 
the d 4 or rather to the d 5 + h situation, where the rel- 
evant ti hole state is partly built from the neighboring 



anion wave functions owing to a strong p — d hybridiza- 
tion. 

We also have simulated the Mnc a K-cdge absorption 
spectra in a GaggMniNgg cluster (a 4a x 4a x 3c su- 
percell, corresponding to 1% Mn concentration) focus- 
ing the attention on the Mn electronic configuration: 
3d 4 and 3d 5 . The calculation is conducted within the 
multiple-scattering approach implemented in fdmnes 69 
using muffin-tin potentials, the Hcdin-Lunqvist approxi- 
mation for their energy-dependent part, a self-consistent 
potential calculation for enhancing the accuracy in the 
determination of the Fermi energy and the in-plane polar- 
ization (£ 1 c). Despite it is common practice to report 
convoluted spectra to mimic the experimental resolution, 
we find out that this procedure can arbitrary change the 
layout of the pre-edge peaks and for this reason it is pre- 
ferred to show non-convoluted data [Fig. 8(6))]. Regard- 
ing the fine structure of the simulated spectra, we have 
a good agreement with experimental data, confirming 
the Mnca incorporation as found by the EXAFS anal- 
ysis (Sec. IV C). On the other hand, the simulated pre- 
edge features need a further investigation: the experi- 
mental intensity of Ai and A2 and the position of A3 are 
not properly reproduced. This could be due to some ne- 
glected effects in the employed formalism, as explained 
in Rcf. 61, where the two peak structure was reproduced 
theoretically within a more elaborated model. 



VI. MAGNETIC PROPERTIES 

A. SQUID results 

We investigate both the temperature dependence of 
the magnetization M at a constant field M(T) and the 
sample response to the variation of the external field at 
a constant temperature M(H). The same experimen- 
tal routine is repeated for both in-plane and out-of-plane 
configurations, that is with magnetic field applied per- 
pendicular and parallel to the hexagonal c-axis, respec- 
tively. In Fig. 9 representative low temperature M (H) 
data for both orientations arc reported. We note that 
these curves exhibit a paramagnetic behavior with a pro- 
nounced anisotropy with respect to the c-axis of the crys- 
tal. This indicates a nonspherical Mn ion configuration, 
expected for aL^O state. At the same time we report an 
absence of any ferromagnetic-like features, that — on the 
other hand — are typical for (Ga,Fe)N layers 28 at these 
concentrations of the magnetic ions, supporting the ab- 
sence of crystallographic phase separation in our layers, 
as suggested by the SXRD and HRTEM studies. The 
same finding additionally indicates that both chemical 
phase separation (spinodal decomposition) and medium- 
to-long range ferromagnetic spin-spin coupling are also 
absent in this dilute layers. The latter allows us to treat 
the Mn ions as completely non-interacting, at least in 
the first approximation. The solid lines in Figs. 9 and 
10 represent fits to our experiential data on the param- 
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agnctic response of non-interacting Mn + ions (L = 2. 
S = 2) with the trigonal crystal field of the wurzite GaN 
structure and the Jahn- Teller distortion taken into ac- 
count (details in Sec. VI B). The overall match validates 
our approach, which, in turn, is consistent with previ- 
ous findings 71 ' 72 that without an intentional codoping, 
or when the stoichiomctry of GaN:Mn is maintained, Mn 
is occupying only the neutral Mn 3+ acceptor state. Inter- 
estingly, all theoretical lines in Figs. 9 and 10 are calcu- 
lated employing only one set of crystal field parameters 
(as listed in Table III) having the Mn 3+ concentration 
n.Mn 3 + as the only adjustable parameter for each indi- 
vidual layer. In Fig. 11 the n Mn !+ values as a function 
of the manganese precursor flow rate are given together 
with the total Mn content iMn as determined by SIMS. 

However, there are hints that the interaction between 
Mn spins may play a role for x > 0.6%. In Fig. 12 
the M(H) normalized at high field (H = 50 kOe) to 
their in-plane values, arc reported. The fact that the 
shape of their magnetization curves is independent of x 
for x < 0.6% means that the interactions between Mn 
ions are unimportant for these dilutions. On the other 
hand, the M (H) for a layer with x = 0.9% (490A) secedes 
markedly from the curves for samples with x < 0.6%, in- 
dicating that supposedly ferromagnetic Mn-Mn coupling 
starts to emerge with increasing relative number of Mn 
nearest neighbors in the layers. Nevertheless, due to the 
generally low Mn concentration in the considered sam- 
ples, no conclusive statement about the strength of the 
magnetic couplings can be drawn from our magnetization 
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FIG. 9: (Color online) Magnetization measurements at 1.85, 
5, and 15 K of Gai-nMn^N as a function of the magnetic 
field applied parallel (closed circles) and perpendicular (open 
squares) to the GaN wurtzite c-axis. The solid lines show 
the magnetization curves calculated according to the group 
theoretical model for non-interacting Mn 3+ ions in wz-GaN. 
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FIG. 10: (Color online) Temperature dependence of the mag- 
netization M for sample 375A (points) at H — 10 kOe. The 
solid lines represent the magnetization calculated within the 
group theoretical model of non-interacting Mn 3+ ions in wz- 
GaN. 



data. Interestingly, depending on the very nature of the 
Mn centers both ferromagnetic and/or antiferromagnetic 
d-d interactions can emerge in (Ga,Mn)N. The presence 
of Mn 2+ ions essentially leads to antiferromagnetic su- 
perexchange, as in II-Mn-VI DMS, where independently 
of the electrical doping, the position of the Mn d-band 
guarantees its 3d 5 configuration. Significantly, the same 
antiferromagnetic d-d ordering and paramagnetic behav- 
ior typical for S = 5/2 of Mn 2+ was reported in n-type 
bulk (Ga,Mn)N samples containing as much as 9% of 
Mn. 73 On the other hand, calculations for Mn 3+ within 
the DFT point to ferromagnetic coupling 74,75 and, exper- 
imentally, a Curie temperature Tq ~ 8 K was observed 
in single-phase Gai-^Mn^N with x ~ 6% and the ma- 
jority of Mn atoms in the Mn 3+ charge state. 11,76 Our 
experimental data seems to support these findings and 
to extend their validity towards the very diluted limit. 
Finally, we remark that the carrier-mediated ferromag- 
netism can be excluded at this stage due to the insulating 
character of the samples, confirmed by room temperature 
four probe resistance measurements and consistent with 
the mid-gap location of the Mn acceptor level. 

The observations presented here point to an uniqueness 
of Mn in GaN. The fact that Ga^^Mn^N with x < 1% 
is paramagnetic without even nanometer-scale ordering 
should be contrasted with GaN doped with other TM 
ions. Depending on the growth conditions, the TM solu- 
bility limit is rather low and typically, except for Mn, it 
is difficult to introduce more than 1% of magnetic im- 
purities into randomly distributed substitutional sites. 
For example, the solubility limit of Fe in GaN has been 
shown to be i k 0.4% at optimized growth conditions 
(see Ref. 28), but signatures of a nanoscale ferromag- 
netic coupling are observed basically for any dilution. 28 
The relatively large solubility limit of Mn in GaN, in turn, 
has a remarkable significance in the search for long-range 
coupling mediated by itinerant carriers. 1,6 Not only it 



11 



„ 5 

CO 

'E 

S° 4 
g 

J= 2 
c 

(D 
O 

§ 1 

O 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — l~~ 



SIMS 
♦ SQUID 



i 



• i 



B i*t 



j i i ■ X ' i i i i i i i i i i i i i i— 



- 1.0 



- 0.5 >< 



100 200 300 400 500 
MeMnCp [ seem ] 



0.0 



FIG. 11: (Color online) Mn concentration n-Mn obtained from 
magnetization measurements (circles - series A, diamonds - 
series B) and SIMS (squares - series A) as a function of the 
Mn precursor flow rate. 
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FIG. 12: (Color online) Magnetization curves at T=1.85K of 
five Gai-zMn^N samples normalized with respect to their in 
plane magnetization at H — 50 kOe. 



lets foresee a high concentration of substitutional Mn — 
important for the long-range ordering — but it can ensure 
that the effects brought about by carriers are not masked 
by signals from nanocrystals with different phases. 



B. Magnetism of Mn ions - theory 

The Mn concentration in our samples is then x < 1% 
as evaluated by means of various characterization tech- 
niques (subsection VA), implying that most of the Mn 
ions (> 90%) have no nearest magnetic neighbors. There- 
fore, the model that considers the Mn ions as single, non- 
interacting magnetic centers should provide a reasonable 
picture. To describe the Mn 3+ ion we follow the group 
theoretical model developed for Cr 2+ ion by Vailin 77,78 
and then successfully used for a Mn-doped hexagonal 
GaN semiconductor. 79,80 It should be pointed out that 
symmetry considerations cannot discriminate between 
the d 5 +hole and d A many-electron configurations of the 
Mn ions, therefore the presented model should be ap- 



plicable to both configurations. Through this section, 
the capital letters T, (i = 1, 2), E denote the irreducible 
representations of the point group for the multielectron 
configurations in contrast to the single electron states in- 
dicated by the small letters d, e, t%. 

We consider a Mn ion that in its free state is in the 
electronic configuration d 5 s 2 of the outer shells. When 
substituting for the group III (s 2 p) cation site, Mn gives 
three of its electrons to the crystal bond and assumes 
the Mn 3+ configuration. In a tctrahcdral crystal field, 
the relevant levels are five-fold degenerate with respect 
to the projection of the orbital momentum and are splet 
by this field and by hybridization with the host orbitals 
into two sublevels e and t% with different energies. In the 
tctrahcdral case the e states lie lower than the t-x states. 
This fact can be understood by analyzing the electron 
density distribution of the t xy , t% z , t% x and e x ~ v , e z 
levels. The density of the state extends along the di- 
rection toward the N ligand anions, while the e orbital 
has a larger amplitude in the direction maximizing the 
distance to the N ion and due to the negative charge 
of the N anions, the ti energy increases. However, the 
relevant, i.e., the uppermost ti state may actually orig- 
inate from orbitals of neighboring anions, pull out from 
the valence band by the p-d hybridization. 81 If the sys- 
tem has several localized electrons, they successfully oc- 
cupy the levels from the bottom, according to Hund's 
first rule, and keep their spins parallel. By considering 
the full orbital and spin moments, the Mn 3+ center can 
be described through the following set of quantum num- 
bers (LrriLSms) with L = 2 and S = 2. However, we 
underline that this procedure can be used only if the 
intra-atomic exchange A ex interaction is larger than the 
splitting between the e and £2 states Ac7_f = E tr , — E e 
(A ex > Acf)- After this, the effect of the host crystal 
is taken into account as a perturbation like in the single 
electron problem. One forms first 2L + 1 wave functions 
for the n-electron system determined by the Hund's rule, 
calculates the matrix elements for these states and deter- 
mines the energy level structure. In this way, the impu- 
rity ions states are found and classified according to the 
irreducible representations of the crystal point group and 
characterized by the set (TM Sms) of quantum numbers, 
with M the number of the line of an irreducible represen- 
tation r = A\, Ai, E, Xi, T2 of the corresponding point 
group. In the case of a Mn 3+ (L = 2, S = 2) ion in a 
tctrahcdral environment the ground state corresponds to 
the 5 T2(e 2 t 2 ) configuration with two electrons in the e 
and two electrons in the level. The ground state is 
three-fold degenerate, since there are three possibilities 
to choose two orbitals from three ti orbitals. The first 
excited state for the Mn 3+ ion is 5 E{e 1 tl) (see Ref. 82). 

The energy structure of a single ion in Mn 3+ charge 
state can be described by the Hamiltonian 



H — Hcf + Hjt + Htr + Etso + #b, (1) 
where H CF = -2/3S 4 (O2-20\/26|) gives the effect of 
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TABLE III: Parameters of the group theoretical model used 
to calculate the magnetization of Gai-^Mn^N. All values are 
in meV. 
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B° 2 


Bl 
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11.44 


4.2 


-0.56 


-5.1 


-1.02 


5.0 


10.0 





a host having tetrahedral Td symmetry, Hjt = B 2 Q\ + 
B®Q\ is the static Jahn- Teller distortion of the tetragonal 
symmetry, Htr = -B O° + B\0\ represents the trigonal 
distortion along the GaN hexagonal c-axis, that lowers 
the symmetry to C^y, -ffso = XLS corresponds to the 
spin-orbit interaction and Hb = /is(L + 2S')B is the Zee- 
man term describing the effect of an external magnetic 
field. Here 0, O are Stevens equivalent operators for a 
tetragonal distortion along one of the cubic axes [100] 
and trigonal axis [111] || c (in a hexagonal lattice) and 

B^, Xtt, and Xte are parameters of the group theo- 
retical model. As starting values we have used the pa- 
rameters reported for Mn 3+ in GaN:Mn,Mg 80 which de- 
scribe well the magneto-optical data on the intra-center 
absorption related to the neutral Mn acceptor in GaN. 
Remarkably, only a noticeable modification (about 10%) 
of Xtt and B\ has been necessary in order to reproduce 
our magnetic data (the remaining parameters are within 
3% of their previously determined values.) Actually, the 
model with the parameter values collected in Table III de- 
scribes both the magnetization M(H) and its crystalline 
anisotropy (Figs. 9 and 10) as well as the position and 
the field-induced splitting of optical lines. 80 

The ground state of the Mn 3+ ion is an orbital and 
spin quintet 5 D with L = 2 and S = 2. The term Hcf 
splits the 5 D ground state into two terms of symmetry 
5 E and 5 T 2 (ground term). The 5 E — 5 T 2 splitting is 
&cf = I2OB4. The nonspherical Mn 3+ ion undergoes 
further Jahn- Teller distortion, that lowers the local sym- 
metry and splits the ground term 5 T 2 into an orbital sin- 
glet 5 B and an higher located orbital doublet 5 E. The 
trigonal field splits the 5 E term into two orbital singlets 
and slightly decreases the energy of the 5 B orbital sin- 
glet. The spin-orbital term yields further splitting of the 
spin orbitals. Finally, an external magnetic field lifts all 
of the remaining degeneracies. 

For the crystal under consideration, there are three 
Jahn- Teller directions: [100], [010] and [001] (center A, 
B, C respectively). 79,80 It should be pointed out that the 
magnetic anisotropy of the Mn 3+ system originates from 
different distributions of nonequivalent Jahn- Teller cen- 
ters in the two orientations of the magnetic field and the 
hexagonal axial field Htr along the c-axis. This picture 
of Mn in GaN emphasizing the importance of the Jahn- 
Tellcr effect, which lowers the local symmetry and splits 
the ground term T 2 into an orbital singlet and a doublet, 
is in agreement with a recent ab initio study employing 
a hybrid exchange potential. 59 

The energy level scheme of the Mn 3+ ion is calculated 



through a numerical diagonalization of the full 25 x 25 
Hamiltonian (1) matrix. The average magnetic moment 
of the Mn ion m = L+2S (in units of /j,b) can be obtained 
according to the formula: 

< m >= Z~ X {Z A < m > A +Z B <m> B +Z C < m > c ), 

(2) 

with Zi (i — A, B or C) being the partition function 
of the i-th center, Z = Za + Zb + Zq and 

i £f=i <<Pi\L + 2S\tp 3 > cxp(-E*/k B T) 

< 111 > = — M I 

Ef =1 exp(-^/fc B T) 

. (3) 

where E % - and ipj are the j-th energy level and the eigcn- 
state of the Mn 3+ ion i-th center, respectively. As already 
mentioned, the Mn concentration in our samples is rel- 
atively small x < 1%. Therefore, the model assuming a 
system of single Mn ions provides a reasonable descrip- 
tion of the magnetic behavior. The macroscopic magneti- 
zation M, shown in Figs. 9 and 10, can then be expressed 
in the form 

M = hb < m > n Mn , (4) 

where nMn = Nmii/V is the Mn concentration and 
A^Mn the total number of Mn ions in a volume V. 



C. Search for hole-mediated ferromagnetism 

As already mentioned, according to the theoretical pre- 
dictions within the p — d Zener model, 1,6 RT ferromag- 
netism is expected in single-phase (Ga,Mn)N and related 
compounds, provided that a sufficiently high concentra- 
tion of both substitutional magnetic impurities (near 5% 
or above) and valence-band holes will be realized. The 
latter condition is a more severe one, as the high binding 
energy of Mn acceptors in the strong coupling limit leads 
to hole localization. 

Surprisingly, RT ferromagnetism in p-type 
Gai^MrijN with a Mn content as low as 
x w 0.25% was recently reported. 13 The investi- 
gated modulation-doped structure consisted of a 
(Ga,Mn)N/(Al,Ga)N:Mg/GaN:Si (i-p-n) multilayer 
and a correlation between the ferromagnetism of the 
(Ga,Mn)N film at 300 K and the concentration of 
holes accumulated at the (Ga,Mn)N/(Al,Ga)N:Mg 
interface was shown. The interfacial hole density was 
controlled by an external gate voltage applied across 
the p-n junction of the structure, and a suppression 
of the FM features — already existing without the gate 
bias — took place for a moderate gate voltage applied. 
Apart from a high value of Tq, a puzzling aspect of 
the experimental results is the large magnitude of the 
spontaneous magnetization, 75 ucmu/cm 2 . 13 Since the 
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holes are expected to be accumulated in a region with 
a thickness of the order of 1 nm, the reported magnetic 
moment is about two orders of magnitude larger than 
the one expected for ferromagnetism originating from an 
interfacial region in (Ga,Mn)N with x = 0.25%. 

Nevertheless, we have decided to check the viability 
of this approach that not only seemed to result in high 
temperature FM in GaN:TM, but also allowed the all- 
electrical control of FM. Thus, we have combined the p- 
type doping procedures we previously optimized 27 with 
the growth of the dilute (Ga,Mn)N presented in this work 
to carefully reproduce the corresponding structure. 13 The 
desired architecture of the investigated sample is con- 
firmed by SIMS profiling (see Fig. 13) indicating the 
formation of well defined (Ga,Mn)N/(Al,Ga)N:Mg and 
(Al,Ga)N:Mg/GaN:Si interfaces. However, as shown in 
Fig. 14, no clear evidence of a ferromagnetic-like re- 
sponse is seen within our present experimental resolu- 
tion of « 0.3 /Ltemu/cm 2 . To strengthen the point, we 
note here that the maximum error bar of our results 
(w 0.7 /^emu/cm 2 ) corresponds to about 1/100 of the 
saturation magnetization reported in the assessed exper- 
iment. While the absence of a ferromagnetic response at 
the level of our sensitivity is to be expected, the presence 
of a large ferromagnetic signal found in Ref. 13 in a nom- 
inally identical structure is surprising. Without a care- 
ful structural characterization of the sample studied in 
Ref. 13 by methods similar to those we have employed in 
the case of our layers, the origin of differences in magnetic 
properties between the two structures remains unclear. 



VII. SUMMARY 

In this paper we have investigated Gai-^Mn^N films 
grown by MOVPE with manganese concentration x < 
1%. A set of experimental methods, including SXRD, 
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FIG. 13: (Color online) SIMS depth profiles of our 
(Ga,Mn)N/(Al,Ga)N:Mg/GaN:Si (i-p-n) structure. 
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FIG. 14: (Color online) Room temperature magnetic sig- 
nal from the (Ga,Mn)N/(Al,Ga)N:Mg/GaN:Si structure. For 
completeness, results of both in-plane and out-of-plane orien- 
tations are shown. Diamagnetic and paramagnetic contribu- 
tions have been compensated. 



HRTEM, and EXAFS, has been employed to deter- 
mine the structural properties of the studied material. 
These measurements reveal the absence of crystallo- 
graphic phase separation and a Ga-substitutional posi- 
tion of Mn in GaN. The findings demonstrate that the 
solubility of Mn in GaN is much greater than the one of 
Cr (Ref. 83) and Fe (Ref. 18) in GaN grown under the 
same conditions. Nevertheless, for the attained Mn con- 
centrations and owing to the absence of band carriers, the 
Mn spins remain uncoupled. Accordingly, pertinent mag- 
netic properties as a function of temperature, magnetic 
field and its orientation with respect to the c-axis of the 
wurtzite structure can be adequately described by the 
paramagnetic theory of an ensemble of non-interacting 
Mn ions in the relevant crystal field. Our SQUID and 
XANES results point to the 3+ configuration of Mn in 
GaN. However, the collected information can not tell be- 
tween d 4 and d 5 + h models of the Mn 3+ state, that is 
on the degree of hole localization on the Mn ions. A 
negligible contribution of Mn in the 2+ charge state in- 
dicates a low concentration of residual donors in the in- 
vestigated films. Our studies on modulation doped p- 
type Gai_j;Mn x N/(Ga,Al)N:Mg heterostructures do not 
reproduce the high temperature robust ferromagnetism 
reported recently for this system. 13 
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